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ABSTRACT 

We report observations of the reactive molecular ions OH + , ILO*, and HjO + towards Orion KL with Herschel/HIFI. All three N = 1 - fine- 
structure transitions of OH + at 909, 971, and 1033 GHz and both fine-structure components of the doublet ortho-H 2 + l n - m transition at 
1115 and 1139 GHz were detected; an upper limit was obtained for H30 + . OH + and H20 + are observed purely in absorption, showing a narrow 
component at the source velocity of 9 km s -1 , and a broad blueshifted absorption similar to that reported recently for HF and para-K^O, and 
attributed to the low velocity outflow of Orion KL. We estimate column densities of OH + and H20 + for the 9 km s -1 component of 9 d 



, and those in the outflow of 1.9 + 0.7 x 10 13 cirr 2 and 1.0±0.3x 10 13 cnr 2 . Upper limits of 2.4x 10' 



: and8.7xl0 12 cm 



and 7: 

were derived for the column densities of ortho and para-H^O* from transitions near 985 and 1657 GHz. The column densities of the three ions 
are up to an order of magnitude lower than those obtained from recent observations of W3 1C and W49N. The comparatively low column densities 
may be explained by a higher gas density despite the assumption of a very high ionization rate. 

Key words. ISM: abundances — ISM: molecules 



1. Introduction 

The Heterodyne Instrument for Far Infrared (HIFI) on the 
Herschel Space Observatory^ provides a unique opportunity to 
fully assess the first steps of the oxygen chemistry in a wide 
variety of sources. Initial HIFI observations quickly detected 
widespread absorption by OH + and H20 + toward the star form- 
ing regions DR21, W31C, and W49N (Ossenkopf et al. 2010; 
Gerin et al. 2010; Neufeld et al. 2010, this issue). Prior to the 
HIFI observations, OH + had only been detected in absorption 
toward Sgr B2(M) (Wyrowski et al. 2010). Similarly, previous 
observations of H20 + were limited to its detection in comet tails 
(e.g., Herzberg & Lew 1974; Wehinger et al. 1974), demonstrat- 
ing the importance of photoionization in producing this ion in the 
absence of H2. And until recently, only upper limits had been re- 
ported on the column density of H20 + in the diffuse interstellar 
gas (Smith, Schempp, & Federman 1984). 

By contrast, the recent HIFI detections of OH + and H20 + in 
warm diffuse gas with a fairly small fraction of molecular hydro- 
gen, elucidated the role of + in initiating the oxygen-hydrogen 
chemistry. This chemistry is thought to begin with the produc- 
tion of H + and H3 via cosmic ray or X-ray ionization of hydro- 



1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



gen, followed by charge transfer to produce + . Rapid hydrogen 
abstraction reactions of + with H2 then yield OH + and H20 + , 
and terminate with the production of H30 + . In diffuse molecu- 
lar clouds, which have high electron abundances, the H30 + is 
destroyed via dissociative recombination to yield OH and H2O. 
In dense molecular clouds, both the ionization fraction and the 
atomic hydrogen abundance are comparatively lower, and the 
sequence of reactions, expected to start at Hi and OH + , yields a 
larger abundance of H30 + . This picture is probably overly sim- 
plistic for molecular clouds such as Orion KL, which are com- 
posed of both diffuse and dense gas. 

Orion KL is the brightest infrared region in the Orion- 
Monoceros molecular cloud complex located less than 500 pc 
from the sun (Menten et al. 2007). In the foreground of Orion KL 
is the Orion Nebula, an HII region known to contain a cluster 
of thousands of young stars which produce a substantial flux of 
X-ray photons (Getman et al. 2005). Molecular line studies re- 
veal three main regions in Orion KL: i. a core of very dense 
and hot gas (n ~ 10 7 cm -2 , T ~ 200 K); ii. cool, quiescent 
gas between systemic velocities of 8 km s _1 and 10 km s , 
surrounded by high-velocity outflows (> 100 km s _1 ); and iii. 
a highly inhomogeneous and turbulent outflow source contain- 
ing both high-velocity (> 30 km s ) and low-velocity (~ 18 
km s" 1 ) gas (Blake et al. 1987; Genzel & Stutzki 1989; O'Dell 
et al. 2008). 
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In this Letter we report the detection of absorption lines of 
OH + and H20 + , and an upper limit on the column density of 
H30 + toward Orion KL. In addition to molecular absorption at 
a systemic velocity of 9 km s _1 , these observations find broad 
blueshifted absorption by OH + and H20 + extending to large neg- 
ative velocities. This is consistent with previously observed lines 
of H 2 with ISO (Lerate et al. 2006), as well as those of HF 
and para-H^O detected recently with HIFI, and attributed to 
the low-velocity molecular outflow (Phillips et al. 2010). 



2. Observations and data reduction 

The observations were done in March 2010 as part of the Key 
Program Herschel/HIFI Observations of Extraordinary sources: 
The Orion and Sagittarius Star-forming Regions (HEXOS). 
The dual beam switch (DBS) observing mode was used, with 
the DBS reference beams lying approximately 3' east and 
west of the Orion KL position a/2000 = 5''35" ! 14.3 s and 
<5y2000 = -5 °22'33.7". Spectra were taken with the Wide Band 
Spectrometer (WBS) with a Nyquist-limited frequency resolu- 
tion of approximately 1 . 1 MHz over a 4 GHz wide IF band; the 
HIFI beams in bands 4, 5, and 6 have half -power beam widths of 
21", 19", and 13" and main beam efficiencies of 0.670, 0.662, 
and 0.645 (HIFI Observers' Manual, v 2.0). The spectra were 
reduced through the standard Herschel Pipeline to Level 2 using 
HIPE version 2.4 (Ott 2010). The double sideband (DSB) spec- 
tra so obtained were then deconvolved (Comito & Schilke 2002) 
to single sideband (SSB) spectra using the doDeconvolution task 
in HIPE. The SSB spectra were converted to the FITS format and 
analyzed with the CLASS90 package. Although two orthogonal 
polarizations were observed simultaneously, only spectra from 
the H polarization in bands 4a and 6b and the V polarization 
in band 5a are shown, because of the smaller standing waves in 
these polarizations. 



Table 1. Spectroscopic parameters of the observed transitions. 
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Notes. {a) Dipole moments (p): 2.256 D (OIL; Werner, Rosmus, 
& Reinsch 1983); 2.37 D (H 2 + ; Wu et al. 2004); 1.44 D (H 3 + ; 
Botschwina, Rosmus, & Reinsch 1985). Frequencies from: (6) Miiller 
et al. (2005); (r) Miirtz et al. (1998); W) Yu et al. (2009) 
(e) / T A dv < 0.482 K km s" 1 for the 984.7 GHz line, and < 2.412 
K km s" 1 for the 1657.2 GHz line. 

(/) Blended with a strong 2 12 - l u ort/jo-Hf O line at 1655831 MHz. 



3. Spectroscopy 

The spectroscopy of OH + , H20 + , and H30 + has been discussed 
in detail in the recent detection papers (Ossenkopf et al. 2010; 
Gerin et al. 2010). Here we summarize the essential aspects of 
the rotational spectra of these ions. The OH + ion has a 3 2~T elec- 
tronic ground state, the two unpaired electron spins (5 = 1) 
yielding three components of the N = 1 - transition. The nu- 
clear spin of the hydrogen atom (Jh = 1/2) further splits each 
component into hyperfine components. The H 2 + ion has Ci v 
symmetry and a 2 B\ ground state which results in the lowest 
level having ortho symmetry. The spin of the unpaired electron 
(S = 1/2) results in two fine-structure components, each ex- 
hibiting a complex hyperfine pattern due to the spins of the two 
equivalent hydrogen nuclei (Jh = 1/2). Rotational spectroscopy 
of H20 + is limited to two laser magnetic resonance (LMR) stud- 
ies (Strahan et al. 1986; Miirtz et al. 1998). Here, we adopt the 
values of Miirtz et al., which we and others have checked inde- 
pendently to be accurate to about 2 km s _1 in equivalent radial 
velocity (see Neufeld et al. 2010 and Schilke et al. 2010, this 
volume). H30 + is a closed-shell symmetric top molecule with a 
large amplitude inversion near 1 .65 THz, resulting in a spectrum 
similar to NH3 with transitions between symmetric and antisym- 
metric inversion states (Yu et al. 2009). Table 1 lists the observed 
transitions of the three ions, along with their line strengths and 
spontaneous emission rates. 



4. Results 

Figure Q] shows the absorption lines of OH + and H20 + toward 
Orion KL, as well as lines of HF and para-H^O for comparison. 
The strongest hyperfine components of OH + and H20 + appear 
at the source velocity of 9 km s _1 , which matches well that of the 
HF line in Orion KL. Additionally, lines of both ions show broad 
blue absorption wings extending to about -75 km s _1 , more ex- 
tended than the HF absorption, but comparable to that of para- 

Hj 8 ( 80 km s~'). We attribute the extended absorption of 

the ions to originate mainly from the low velocity molecular out- 
flow. We failed to detect any emission or absorption from H30 + , 
and discuss the non-detection in § [5] 

The high density of molecular lines in Orion KL makes 
contamination by unrelated lines a common problem. The ab- 
sorption lines detected here are blended with weak to mod- 
erately strong emission lines of abundant "weeds", including 
CH3OH and SO2. Efforts are underway to model and remove 
the emission from the contaminants by a method similar to that 
of Phillips et al. (2010); in the interim, the following approach 
was taken. 

To better gauge the absorption, the contaminants were 
masked and intensities interpolated across the masked channels 
(Fig- CJ. The velocity-integrated optical depths of the ionic lines 
were obtained by normalizing the SSB spectra with the con- 
tinuum and integrating over the velocity ranges for the source 
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of 1.9 



: 0.7 x 10 13 cm" 2 and 1.0 + 0.3 x 10 13 cm" 2 . The column 
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Fig. 1. Lines of OH + , ortho-H 2 + , and H 3 + in Orion KL, com- 
pared with those of HF (J = 1-0) and para-H^O. The dashed 
vertical red line is at the systemic velocity of 9 km s , and the 
solid horizontal red lines indicate the continuum level in each 
spectrum. Solid green lines indicate the channels over which the 
interpolation was done. The HF spectrum, adapted from Phillips 
et al (2010), shows a broad absorption (black histogram) after 
modeling and removal of contaminating lines of CH3OH and 
SO2 (blue histogram). The contaminants in the 971.8 GHz spec- 
trum are CH3OH (v, = 0) 26j -25j (A) and (y, = l)22 6 -22 5 (E) 
at -28 km s _1 and -40 km s _I ; that in the 1 1 15.2 GHz spectrum 
is S0 2 19 5; i5 - 19 2 ,is at -55 km s -1 . The 1657.2 GHz spectrum 
shows absorption by CH (N — 2 - 1) at 60 km s . 



and the outflow, the interpolation yielding errors of 20% - 30%. 
On the assumptions that the absorption covers the source com- 
pletely, and the molecules are in the lower state, the total column 
density (N) was then derived using the expression: 



rdv (km s 1 ) 



N, 



(1) 



where A„/ is the spontaneous emission rate, g u and gi are the 
upper and lower state degeneracies, and A is the transition wave- 
length. 

We estimate column densities of OH + and H20 + at 9 km s _1 
of 9 + 3 x 1 1 2 cirT 2 and 7 + 2 x 1 12 cirT 2 , and those in the outflow 



densities of OH + are more than an order of magnitude lower, 
and those of H20 + are 2 - 6 times lower than toward W3 1C and 
W49N (Gerin et al. 2010; Neufeld et al. 2010, this issue). From 
the least congested spectra of H 3 + at 984.7 and 1657.2 GHz 
(see Table 1), and an assumed excitation temperature of 100 K, 
we derive 3<x upper limits of 2.4 x 10 12 crrT 2 and 8.7 x 10 12 cm -2 
for the column density of ortho and para-Hj,0 + , nearly an order 
of magnitude lower than in W31C (Gerin et al. 2010). 

The abundance ratios of the three ions in Orion KL can be 
compared to the same ratios observed in W31C and W49N. The 
OH + /H 2 + ratio is found to be 1.3 + 0.6 in the source and 1.8 ± 
0.8 in the outflow. This ratio is 2 - 15 times lower than that 
measured toward W31C and W49N. The lower limit of 1.4 for 
the H20 + /H 3 + ratio, however, is nearly 2 times larger than in 
W31C. 



5. Discussion 

The column densities of OH + , H20 + , and H 3 + in Orion KL 
differ markedly from those in the diffuse gas toward W31C and 
W49N. In contrast with W31C and W49N, OH + and H 2 + are 
significantly more abundant relative to H 3 + , for which we are 
only able to obtain an upper limit. The absolute column densi- 
ties of OH + and H20 + are also lower compared with W31C and 
W49N. A likely explanation for the low column densities of the 
three ions is that they are present in fairly dense material, both in 
the quiescent gas and the outflow. Unlike the quiescent gas, the 
Orion KL outflow is exposed to a strong ionizing flux from the 
foreground HII region; the enhanced ionization flux enhances 
the formation of ions, but the resultant large fractional ioniza- 
tion leads to a fast and efficient removal of molecular ions by 
dissociative recombination with electrons. 

The observed velocity profiles of OH + and H20 + in 
Orion KL support the above conclusion. As Fig. Q] shows, the 
OH + and H20 + absorption tracks the HF absorption to veloci- 
ties of about -45 km s _1 . This absorption also seems to follow 
closely, to about -80 km s _1 , the para-H^O absorption in the 
outflow, suggesting that like HF and para-H^O, OH + and H20 + 
probably exist mainly in the low velocity outflow (Phillips et al. 
2010). In fact, the molecular outflow accounts for over half of 
the observed column density of OH + and H20 + . 

The conditions required to explain our observations may be 
more extreme than one might suppose. First, the molecular ions 
probably reside in gas of lower density (n < 10 s cirT 3 ) than that 
necessary to thermally excite the observed transitions — these 
have high spontaneous emission rates (> 10~ 2 s ; Table 1), and 
hence large critical densities (10 7 — 10 9 cm" 3 ) . This is supported 
by the observation that OH + and H20 + are seen only in absorp- 
tion. Second, the temperatures in the outflow gas are probably 
high. 



We consider two scenarios in which the ions may be formed 
in the low velocity outflow. In the first, a large radiation flux im- 
pinges directly on the Orion KL outflow, which contains large 
water abundances (Melnick et al. 2010). The far UV flux that il- 
luminates this gas can have values approaching 4 x 10 4 times 
the average interstellar radiation field (Walmsley et al. 2000; 
Young Owl et al. 2000). In addition, the central region of the 
Orion Nebula has numerous sources of energetic X-ray photons 
(Getman et al. 2005; Preibisch et al. 2005), which can contribute 
to the surface ionization of this photon-dominated region (PDR). 
We estimate that at Ay = 1 into the PDR, the ionization rate 
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Under these conditions, water can undergo 
photoionization to form H20 + directly, enhancing the abundance 
of this species. 

In the second scenario, the outflow penetrates the extended 
foreground HII region. The abundant H + can now undergo 
charge exchange with H2O to yield HaO + . In either scenario, the 
high electron density probably results in a net reduction in the 
abundances of molecular ions, consistent with the observations: 
low column densities of OH + and H20 + , and the upper limit for 
H 3 + . 

We have attempted to model the first scenario using the 
Meudon PDR code (Le Petit et al. 2006). However, the model 
suffers difficulties while reproducing the observed column den- 
sities of the three ions. First, it requires a relatively low gas 
density (n ~ 10 3 cm 3 ) in regions where OH + and H20 + are 
produced, as larger assumed densities yield too much H30 + . 
Second, it requires a very large ionization rate (£ > 1 - 2 x 10~ 14 
s _1 ) to maintain a ratio of atomic to molecular hydrogen near 
unity; otherwise, too much H?0 + is once again produced. The 
two parameters are nearly an order of magnitude different from 
others inferred from previous observations: n > 10 4 cm 4 and 
£ < 10~ 14 s" 1 (Genzel & Stutzki 1989; Lerate et al. 2008; 
Muench et al. 2008 and references therein). Nevertheless, a re- 
cent study on molecular hydrogen rotational excitation in the 
Orion bar infers a cosmic ray ionization rate of 7 x 10~ 14 s _1 
(Shaw et al. 2009). The same study also invokes warm gas tem- 
peratures of 400-700 K; the lower value is contained in our 
model for the edge of the PDR. A critical evaluation of our 
model awaits further work and a thorough exploration of the pa- 
rameter space, and will be presented in a future paper. 

We are unable to confirm previous tentative detections of 
H 3 + toward Orion KL (Hollis et al. 1986; Wootten et al. 1986; 
Wootten et al. 1991; Phillips et al. 1992; Timmermann et al. 
1996; Lerate et al. 2006). Of these, Phillips et al. (1992) present 
the best evidence: 3 emission lines at 307, 364, and 396 GHz, ly- 
ing 45, 85, and 105 ctrr 1 above ground; but they do not rule out 
the possibility of blends with other lines. The lines we observed 
are at lower energies (see Table 1), and are expected to be as 
strong or stronger than those observed by Phillips et al. (1992). 
The upper limits derived here for ortho andpara-H30 + are more 
than an order of magnitude lower than the column densities re- 
ported by Phillips et al. (1992). Timmermann et al. (1996) re- 
ported detection of the 4^ - 3, line near 70 fim with the Kuiper 
Airborne Observatory, but the velocity of the line differs by more 
than -60 km s from predicted values. Lerate et al. (2006) de- 
tected the 2\ - 1|, 2~ - 1+, and 1~ - 1| lines with ISO: the first, 
near 2.98 THz is 80 km s _1 higher than the systemic velocity of 
9 km s ; the second, near 2.97 THz, is 1 km s higher than 
the frequencies predicted by Yu et al. (2009); and the third, at 
1655835 MHz, covered by our observations, is obscured by a 
strong 2i2 - 1 10 ortho-W^O line at 1655831 MHz. 

6. Conclusions 

Our observations toward Orion KL have found OH + and H20 + 
aborption at the quiescent 9 km s _1 component and extended ab- 
sorption in the low velocity molecular outflow associated with 

2 The surface brightness of the central region (dominated by G [ C) is 
estimated to be 3 x 10 34 ergs s~'pc~ 2 (Feigelson et al. 2005). On the 
assumption that the molecular cloud lies 0.1 pc from this cluster, the 
expressions of Maloney et al. (1996) yield an X-ray ionization rate of 
about 2.8 x lO -16 ^, 1 (where N22 is the hydrogen column density in 
units of 10 22 cnr 2 ). Thus at A v = 1, & ~ 3 X 10~ 15 s" 1 . 



this source. This is, to our knowledge, the first detection of 
these ions toward a source with a large fraction of molecular 
gas. Given the complex and inhomogeneous nature of Orion KL, 
however, there are probably regions where the densities are suf- 
ficiently low and the excitation conditions optimal for these re- 
active ions to exist at detectable levels. Another possibility is 
that depletion of some of the gas-phase species onto the grains 
can result in lower abundances of water, leading to small col- 
umn densities of OH + and H20 + . A surprising observation — 
and one remarkably different from that toward W31C — is the 
non-detection of H30 + . In our model of the outflow, we attribute 
this mainly to a very high ionization rate, which produces an al- 
most equal abundance of atomic and molecular hydrogen at the 
assumed density. 
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